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ABSTRACT 


A program has been written which generates a nozzle contour 
and the complete flow field for two dimensional or axisyrn- 
metric flows designed to exit parallel to the axis at uniform 
pressure. The flow is tha.t of a rotational, non-homentropi c 
gas mixture where viscous effects have been neglected and the 
chemistry is assumed frozen. This report comprises a com- 
plete description of the numerical program developed, the 
analysis being described in Advanced Technology Laboratories 
Technical Report No. 148 entitled, "The Determination of 
Nozzle Contours for Rotational, Non-homentropi c Gas Mixtures." 
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1 • INTRODUCTION ' 

The overall nozzle configuration is depicted in Figure (1). 

The numerical procedure and hence the program, may be divided 
into several basic areas: 

(A) Initial Profile. The data is specified at the 
NPTS data points on the non-characteri s ti c initial line W-|A|v- 

(B) Down-Runni ng Characteristics Emanating from 
Initial Profile. For each data point N on the initial line, a 
C_ characteri s ti c is calculated emanating from the initial line 
and ending on the axis. The first C_ characteristic emanates 
from initial line data point N=2 and the final one from N=NPTS. 

On a C_ characteristics emanating from point N on the initial 
line, there are LMAX grid points, where LMAX=2N-1 . Note, that 
to calculate the C characteristic, emanating from N, only the 
value on the initial line and the previous C characteristic, 
emanating from N-l , are required, hence only two successive C 
lines need be stored. 

(C) Down-Running Characteristics Emanating from 
Prescribed Arc. Since the arc is generally a circle of small 
radius, the grid must be refined in the area of the wall. The 
mesh points are hence subdivided in the vicinity of wall corres- 
ponding to the severity of the turn. In addition, logic is in- 
corporated to force the angle change at the wall between two 
successive down-running characteristics to be less than one degree 
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by further grid subdivision. The calculation along C_ charact- 
eristics proceeds from the wall to the axis, and the procedure 
terminates when the axis pressure falls below design pressure. 
The characteri s ti c WgA^, yielding design pressure on the axis 
is found by interpolation between the two previously calculated 
C_ lines 

(D) Generation of Exit C + Characteristic. Let II de- 
note a data point on the down-runni ng characteri sti c W 2 A 3 and ii 
denote a data point, having the same value of streamf uncti on , 
on the up-runni ng characteri s ti c A 3 W 3 . All properties are known 
at il , while at N, the pressure is design pressure and the flow 
is parallel to the axis. Since both mass fraction (ct.,-) and en- 
tropy are constant along NN, knowing the pressure ratio P^/Pq 
yields the temperature T- by an iterative procedure. The con- 
stancy of stagnation enthalpy then enables the velocity q- to 
be determined. Having properties at the mesh points N, we then 
determine the location of the mesh points by use of the stream- 


function equation. 

Having the exit 

characteristic A^W^ calcu- 

lated, the following information is 

available without having 

to proceed with 

the 

contour design: 



1 . 

nozzle length 



2 . 

area ratio 



3. 

complete profiles 

at exit plane 

(E) 

Contour Calculation. 

The calculation works along 


C_ characteri s ti cs now working from data points on the exit 
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characteristic back to the wail, whose location is unknown. 
The wall point calculation is depicted in Figure (2). The 
wall angle 0 C is obtained by an iterative procedure. For 
an assumed value of 0* c , point C is located as the inter- 
section of the streamline from D, with the charac teristi c 
from B. The correct value of 0 C , is that for which the com 
patibility relations applied along both AC and BC yield 
the same pressure P c to within a prescribed tolerance. 

In Section II, the input required is described, while the 
output description is presented in Section III. A descrip- 
tion of the subroutines and functions used in the program 
is presented in Section IV. The specific machine require- 
ments and time estimates are listed in Section V. A de- 
tailed flow chart is presented in Section VI. 




DESCRIPTION OF INPUT 


Col umns 

• Format 

Description 

1 -5 • 

15 

i Type of flow: (0=Two Dimen- 
sional , 1 =Axi symmetri c) 

6-10 

15 

Number of data points on 
initial profile (NPTS) 

11-15 

15 / 

Output indicator (check Des 
cription of Output, Section 
HI) 

1-10 

E10.0 

Ratio of axial coordinate 
of wall at initial station 
to throat radius 

11-20 

E10,0 

Ratio of radial coordinate 
of wall at initial station 
to throat radius 

21-30 

E10.0 

Initial inclination of wall 
; (_in radians) 

31-40 

E10.0. 

Ratio of radius of initial 
circular arc expansion to 
throat radius 

41-50 

E10.0 

Ratio of exit pressure to P 

51-60 

El 0 . 0 

Throat radius 
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Card No. 
3 


Columns 

1-10 

11-20 . 

21-30 


Format 

E10.0 

E 1 0 . 0 

E 1 0. 0 


Description 

Free stream or reference 
Mach number 

Free stream or reference 
temperature (°K) 

Free stream or reference 
average molecular weight 


4 Initial Profile - There are 2 cards required for 

each data point as described be- 
low. Begin inputting data points 
at the axis (Point#!) and proceed 
to the nozzle wall. (Point #NPTS) 


4a 


1-10 E10.0 

11-20 E10.0 


21-30 E10.0. 


31-40 E10.0 


41-50 E10.0 


Ratio of axial coordinate of 
data point to throat radius 
Ratio of radial coordinate 
of data point to throat radi- 
us 

Ratio of velocity at data 
point to free stream velo- 
city 

Ratio of temperature at data 
point to free stream temper- 
a t.ure 

Ratio of pressure at data 
point to free stream pressure 


I 
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Card No. 
4a 


Co 1 Limn s Format 

51-60 El 0.0 


Descr i pti on 

Flow inclination at data 
point (in radians) 


1-10 

E10.0 

Bass 

fraction 

of 

H 

11-20 

E10.0 

Bass 

fraction 

of 

0 

21-30 

E10.0 

Mass 

f racti on 

of 

h 2 o 

31 -40 

E10.0 

Mass 

fraction 

of 

h 2 

41-50 

E10.0 

Mass 

fraction 

of 

02 

51-60 

El 0 . 0 

Mass 

f racti on 

of 

OH 

61-70 . 

E10.0 

Mass 

fraction 

of 
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III . DESCRIPTION OF OUTPUT 

( A ) Output Indicator Equals 1 . If the output indi- 
cator equals 1 the following output is obtained: 

1. Input (control constants, reference condi- 
tions and initial profile) 

2. Down-running characteris ti c line from ini- 
tial wall point to axis ( W i A 2 i n Figure ( 1 } ) 

3. Line number, number of points on the line 
and the wall and axis points of C_ charac- 
teristic lines between the line ( W-| A 2 to 

W g A 3 in Figure ( 1 ) ) 

4. Down-running characteri sti c line with axis * 
point pressure equal to design pressure 

5. Exit characteristic line (beyond which pres- 
sure equals design pressure and flow is par- 
allel to axis) ( w 3 A ^ In Figure (I)) 

6 . Nozzle design points 

7 . Exi t prof i 1 e 

(B) Output Indicator Equals 0 . When the output indi- 
cator is set equal to 0 all the output in paragraph A is ob- 
tained plus a print-out of data along all the C characteristics 
in the flow field, at a selected number of points along them. 
This indicator would only be used when a detailed description 
of the flow field is required such as profiles at intermediary 


axial stations, which could be obtained by an interpolation 
procedure along these C_ characteristics. 
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IV. SUBROUTINES AND FUNCTIONS 


1 . 


(A) Subrouti nes . - 
Name ♦ Descri pti on 

THERM Calculates entropy, S ■ , and enthalpy , 

H.j , from polynomial fits in temperature 


2. ALL Computes subsidiary properties as a 

function of basic variables (Q»T.P*a^) 

C p = 1 C p i a i 

W » ^Zo^/m^y 1 
p * y M 2 W P/T/W ' 

R - R 0 /W 

V - C p /(C p - R/C p J 
M - Q M„(y„R„/y/R/T)^ 
y - tan ^ 1 / ( M 2 - 1 Y i ) 

3. COEFF Sets thermodynamic coefficients as func- 

tions of temperature 


4. THERMO Calculates specific heat, C n , deriva- 

Pi 

tive of specific heat and enthalpy, H^, 
of each species from polynomial fits in 
temperature 
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Name Description 

5. ERROR Outputs program statement number near- 
est error and terminates computer run 

( B ) Functions 

Marne Descri pti on 


1 . 

XM1 

tan (e + 

y) along C + 

characteri Stic 

2. 

XM2 

tan ( e - 

y) along C 

characteri s ti c 

3. 

XM3 

tan (e) 

along streamline 

4. 

FI 

* 

A^ or B-j 

coeff i ci ent 

along C+ charac- 



teri sti c 



5. 

F2 

A g or B* 

coef f i ci ent 

along C + charac- 


teri s ti c 


*No te : P - 0 relationship along C + characteri s t i c : 

A 1 (P C - P A> + e c - e A + A 2 (x 2 ' X A > ' 0 
P - 6 relationship along C_ characteristic: 
B l^ f) C * P B^ “ e C + 0 B + B 2^ X C " X B ^ = 0 
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MACHINE CONTROL CONSIDERATIONS 

( 

1. Machine - program designed for CDC 5600 

2. Loader -PPLOADR 

i 

3. Estimates for run: 

a. Field Length: 

(1) compile - 64000 g locations 

(2) loading - 70000g locations 

b . CP time: variable depending on number of 

points and properties of an initial pro- 
file'and type of flow (approximately 2 
minutes for sample case included) 

c. 10: less than 100 g seconds 

d. Tapes and disks: 

• ' i 

( 1 ) tape 5 - card input 

(2) tape 6 - printed output 

(3) tape 7 - punched output 

(4) no other tapes or disk files are used 

e. Printed output: 7000g lines including list- 

ing of program for output indicator equals *1, 
60000 g lines for output indicator equals 0. 
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FLOW CHART 


Symbols and References in Flow 'Chart 


Basic variables: 


k - 

axial di stance/throat radius 

y 

radial di s tance/throa t radius 

i ■ . 

?■ - 

press ure/free stream pressure 

Q 

vel oci ty/f ree stream velocity 

T 

temperature/f ree stream temperature 

TH 

flow angle 

ALP ( 1 -7 ) 

mass fractions (seven species) 

Subsidiary properties: 

CPX 

specific heat/free stream specific heat 

R 

gas constant (Rq/W) 

GAM ; 

ratio of specific heats 

EM 

Mach number 

XMU 

Mach angle 

RHO 

density /free stream density 

Initial data line - called 

K LINE - example of variable on this 
line is X ( K ) 

New data line - called 

L LINE - example of variable on this 
line is XN(K) 

Streamline 

example of variable on this 


line is XD 

NPTS - number of points on initial profile until 

down-running line giving design pressure at 
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axis, after this, it is number of points 
on up-running characteristic lines 
N - line number from 1 to number of points on 

! • i ni ti al prof i 1 e 


KfyiA X 

- number 

of 

points 

on 

K LINE 

LljiAX 

- number 

of 

po i nts 

on 

L LINE 

DYT 

- one half 

averge 

A y 

on initial profile 

NSAVE 

- number 

of 

points 

on 

both down-running char 


acteristic line giving design pressure at 
axis anu the final up-running characteris- 
tic line 


LF - counter on final up-running characteristic 

line, for values 1 to NSAVE 


T M 1 
Ill i-i-o 
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